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ABSTRACT 

The  direct  approach  of  modeling  the  flow  between  all  blade  pas¬ 
sages  for  each  blade  row  in  the  compressor  is  too  computationally 
intensive  for  practical  design  and  analysis  investigations  with  inlet  dis¬ 
tortion.  Therefore  a  new  simulation  tool  called  the  Turbine  Engine 
Analysis  Compressor  Code  (TEACC)  has  been  developed.  TEACC 
solves  the  compressible,  time-dependent,  3D  Euler  equations  modified 
to  include  turbomachinery  source  terms  which  represent  the  effect  of 
the  blades.  The  source  terms  are  calculated  for  each  blade  row  by  the 
application  of  a  streamline  curvature  code.  TEACC  was  validated 
against  experimental  data  from  the  transonic  NASA  rotor,  Rotor  IB, 
for  a  clean  inlet  and  for  an  inlet  distortion  produced  by  a  90-deg,  one- 
per-revolution  distortion  screen.  TEACC  revealed  that  strong  swirl  pro¬ 
duced  by  the  rotor  caused  the  compressor  to  increase  in  loading  in  the 
direction  of  rotor  rotation  through  the  distorted  region  and  decrease  in 
loading  circumferentially  away  from  the  distorted  region. 

NOMENCLATURE 

Cp  Constant  pressure  specific  heat 

E,  F,  G  Flux  vectors 

e  Internal  energy  and  kinetic  energy  per  mass 

F  Component  of  force 

F  Force 

m  ;  Rate  of  mass  flow 

P  Pressure 

Q  Conservation  variables 

R  Ideal  gas  constant 


Residual  Vector  representing  the  value  of  the  explicit  part  of  the 
implicit  discretized  conservation  equations 

r  Radial  direction 

S  Component  of  the  source  vector 

S  Source  vector 

SW  Rate  of  shaft  work 

T  Temperature 

u  Velocity  in  the  x-direction 

V  Component  of  velocity 

V  Velocity 

V  Volume 

V  Velocity  in  the  y-direction 

w  Velocity  in  the  z-direction 

X,  y,  z  Cartesian  right-handed  coordinate  directions 

Subscripts 

b  Bleed 

t  Total 

INTRODUCTION 

Modem  high-performance  military  aircraft  are  subjected  to  rapid 
flight  maneuvers  which  place  great  operational  demands  on  theft  air- 
breathing  gas  turbine  engines.  One  component  of  the  engine  that  is  par¬ 
ticularly  sensitive  to  the  fluid  dynamic  transients  that  result  from  rapid 
aircraft  maneuvers  is  the  compressor.  The  compressor  should  operate 
in  a  stable  manner  during  all  aspects  of  flight.  However,  rapid  flight 
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transients  cause  the  inlet  to  produce  a  highly  distorted  total 
pressure  flow  field  to  the  compressor  inlet.  High  distortion 
levels  may  cause  the  compressor  to  surge  at  high  rotational 
speeds  or  slip  into  rotating  stall  at  lower  rotational  speeds 
(AIR- 14 19,  1983).  Since  total  pressure  distortion  is  the  pri¬ 
mary  reason  for  reaching  the  engine  stability  limit,  its 
effects  on  system  performance  and  operability  need  to  be 
understood. 

Distortion  imposed  on  a  circumferentially  swirling 
flow  was  shown  by  Greitzer  and  Strand  (1978)  to  have  a 
three-dimensional  (3D)  nature  which  is  fundamental  to  the 
development  of  both  inlets  and  compressors.  Design  or 
analysis  engineers  are  interested  in  understanding  the 
details  of  the  flow  field  to  determine  the  effects  of  inlet 
total  pressure  distortion  on  the  compressor.  One  way  to 
quantify  the  effects  of  distortion  is  to  test  for  that  effect  in  a 
ground  test  facility.  Currently,  the  inlet  and  engine  are  tested  sepa¬ 
rately.  Typically,  the  aircraft  fuselage  is  too  big  to  fit  in  a  wind  tunnel. 
A  forebody  simulator  is  used  in  conjunction  with  the  inlet  to  character¬ 
ize  its  flow  field.  The  forebody  simulator  is  designed  to  produce  a  flow 
field  at  the  inlet  reference  plane  (IRP)  similar  to  the  flow  field  pro¬ 
duced  by  the  aircraft  (Fig.  1). 

Screens  are  constructed  to  capture  the  most  severe  dynamic  pat¬ 
terns  produced  by  the  inlet  and  are  then  placed  in  front  of  the  engine  to 
measure  the  loss  of  stall  margin  produced  by  the  steady-state  inlet  dis¬ 
tortion  (AIR- 1419,  1983).  However,  it  is  expensive  to  instrument  a 
compressor  and  perform  the  necessary  number  of  tests  to  adequately 
understand  the  compressor  flow  field.  Therefore,  numerical  simulations 
have  been  developed  to  support  the  testing  community. 

Parallel  compressor  theory  has  been  used  successfully  to  develop 
an  understanding  of  compressor  performance  with  inlet  distortion.  Par¬ 
allel  compressor  theory  models  the  compressor  with  multiple  circum¬ 
ferential  segments  and  a  uniform  exit  back  boundary  condition.  Parallel 
compressor  theory  is  restricted  to  simple  inlet  distortion  patterns  and  is 
consistently  conservative  at  estimating  the  stability  limit.  Investigators 
have  made  extensive  modifications  to  parallel  compressor  theory 
through  modeling  techniques  to  account  for  the  transfer  of  mass, 
momentum,  and  energy  transfer  between  segments  (Korn,  1974;  Steen- 
ken,  1983;  Mazzaway,  1977;  Longley,  1992;  Kimzey,  1977;  Shahr- 
okhi,  1995).  3D  simulations  have  recently  been  developed  which  auto¬ 
matically  account  for  the  transfer  of  conservation  properties  (Billet, 
1988;  Tan,  1996).  However,  these  investigations  have  been  restricted  to 
incompressible  flow  or  low  speed  compression  systems.  There  is  still  a 
need  to  develop  a  3D  compression  system  to  investigate  the  effects  of 
complex  inlet  distortion  on  a  high  speed  compression  system. 

Therefore,  a  numerical  simulation  of  inlet  distortion  has  been 
developed  and  is  called  the  Turbine  Engine  Analysis  Compressor 
Code  (TEACC).  TEACC  is  a  direct  extension  of  previous  work  by  the 
author  (Hale,  1994)  and  solves  the  compressible  three-dimensional 
(3D)  Euler  equations  over  a  finite-volume  grid  domain  through  each 
blade  row.  The  Euler  equations  are  modified  to  include  turbomachinery 
source  terms  which  model  the  effect  of  the  blades.  The  turbomachinery 
source  terms  are  bleed,  blade  forces  in  the  three  Cartesian  directions, 
and  shaft  work.  The  source  terms  are  calculated  for  each  circumferen¬ 
tial  grid  section  of  each  blade  row  by  the  application  of  a  streamline 
curvature  code.  A  methodology  was  developed  for  distributing  the  tur- 
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Fig.  1.  Engine  screen  developed  to  duplicate  distortion. 

bomachinery  source  terms  axially,  radially,  and  circumferentially 
through  the  bladed  region. 

OVERALL  APPROACH  TO  TEACC 

The  overall  TEACC  development  methodology  is  conceptually 
presented  in  Fig.  2.  A  general  purpose  3D  flow  simulation  computer 
code  (NPARC,  Cooper,  1989)  is  modified  to  accept  turbomachinery 
source  terms  by  semi-actuator  disk  theory.  These  turbomacMnery 
source  terms  are  calculated  using  a  streamline  curvature  code  (SLCC, 
Hearsey,  1970).  TEACC  is  constructed  by  combining  the  technology  of 
solving  the  Euler  equations  (NPARC)  modified  to  include  source  terms 
and  the  technology  of  calculating  the  source  terms  (SLCC)  to  produce  a 
time-dependent  turbomachinery  simulation  with  the  capability  of  ana¬ 
lyzing  inlet  distortion. 

Governing  Equations 

The  governing  equations  used  in  TEACC  are  developed  by  apply¬ 
ing  the  conservation  of  mass,  momentum,  and  energy.  In  turbomachin¬ 
ery  flows,  the  viscous  effects  predominate  mostly  along  the  wall,  mak¬ 
ing  accurate  simulation  of  the  flow  field  away  from  the  wall  possible  by 
using  the  Euler  equations  with  turbomachinery  source  terms  (Longley, 
1992).  The  equations  of  fluid  motion  using  the  thermally  and  calorically 
perfect  ideal  gas  assumption  with  turbomachinery  source  terms  are: 
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d.  Analysis 

i 


b.  SLCC  c.  TEACC 

Fig.  2.  Overall  TEACC  methodology. 

The  volumetric  turbomachinery  source  terms  are:  (1)  bleed  per  vol¬ 
ume,  (2)  forces  in  the  x,  y,  and  z  direction  per  volume,  and  (3)  mte  of 
shaft  wea’k  pei'  volume.  The  technique  for  calculating  these  source 
terms  from  the  SLCC  are  developed  later  in  this  paper. 
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n  =  number  of  nodes.  The  second  criteria  used  to  verify  TEACC’ s  con¬ 
vergence  to  steady-state  is  to  monitor  key  variables  of  interest  and  verify 
they  no  longer  vary  as  TEACC  continues  to  iterate. 

TURBOMACHINERY  SOURCE  TERMS 

Three-dimensional  blade  force  and  rate  of  shaft  work  terms  are 
supplied  by  a  streamline  curvature  code  (SLCC).  The  SLCC  is  based 
on  a  radial  redistribution  of  blade  force  and  shaft  work  producing  an 
“axisymmetric  flow  with  swirl”  in  the  form  of  streamlines.  Necessary 
inputs  include  overall  geometry,  blade  geometry,  and  loss  and  devia¬ 
tion  correlations.  These  correlations  are  obtained  from  experimental 
data.  A  typical  SLCC  grid  illustrated  in  Fig.  3  extends  far  upstream  and 
far  downstream  of  the  compressor. 

Inlet  Boundary 
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TEACC  Boundary  Conditions 

The  time  dependent  Euler  equations  with  source  terms  define  a 
hyperbolic  in  time  and  an  elliptic  in  space  system  of  equations  requiring 
a  full  set  of  boundary  conditions  specified  at  all  boundaries.  The  bound¬ 
ary  conditions  used  with  NPARC  are  explicit.  The  inflow  boundary 
condition  is  based  on  reference  plane  characteristics,  and  the  total  pres¬ 
sure  and  total  temperature  at  the  inlet  are  specified.  Inlet  flow  directions 
are  assumed  to  be  normal  to  the  boundary.  The  exit  boundary  condition 
is  a  variable  static  pressure  capable  of  supporting  the  exit  profile  of 
strong  swirl.  The  exit  static  pressure  is  calculated  by  specifying  a  single 
value  and  imposing  the  static  pressure  profile  of  the  adjacent  upstream 
station  on  the  exit.  The  wall  boundary  conditions  are  assumed  to  be  slip 
waU;  the  normal  velocity  components  are  set  equal  to  zero  at  the  solid 
walls,  A  rotationally  periodic  (wrap-around)  boundary  condition  is  used 
in  the  circumferential  direction,  where  the  circumferential  seam  of  the 
grid  was  overlapped  by  one  circumferential  segment. 

TEACC  Solution  Procedure 

TEACC’ s  and  NPARC ’s  solution  procedures  are  the  same,  solving 
the  governing  equations  with  the  Beam  and  Warrmng  approximate  fac¬ 
torization  algorithm  (Beam,  1976)  to  obtain  a  flow-field  solution.  The 
above  equations  are  consistently  discretized  with  central  differences  to 
produce  an  implicit  algorithm.  The  governing  equations  were  cast  in 
strong  conservation  form  to  keep  global  conservation,  even  if  strong 
shocks  are  present.  The  technique  of  Pulliam  and  Steger  (1980)  to 
implement  the  Beam  and  Warming  implicit  algorithm  was  incorporated 
to  diagonalize  the  implicit  matrices,  requiring  a  penta-diagonal  solver. 
Since  this  is  a  central  difference  algorithm,  artificial  dissipation  is 
included  for  stability.  NPARC  uses  the  Jameson-style  (1981)  improved 
shock-capturing  artificial  dissipation.  TEACC,  as  well  as  NPARC,  uses 
the  time  stepping  method  to  achieve  steady-state.  There  are  two  criteria 
used  to  verify  TEACC’s  convergence  to  steady-state.  First,  the  L2 
norm  of  the  residual  is  dropped  to  machine  zero  defined  as: 

iResidualL  =  /T  V  Residual^^; 


Fig.  3.  Original  SLCC  with  overall  geometry. 

Axial  stations  are  chosen  to  include  the  inlet  and  exit  of  flie  blades 
and  other  convenient  stations  to  allow  a  proper  resolution  of  the  flow 
field.  Boundary  conditions  are  required  around  the  computational 
domain,  with  the  outer  casing  and  inner  center-body  treated  as  stream¬ 
lines.  The  boundary  conditions  are  inlet  total  pressure,  inlet  total  tem¬ 
perature,  and  Inlet  swirl  angle  with  curvature  set  equal  to  zero  at  both 
the  inlet  and  exit.  A  particular  operating  point  is  established  by  specify¬ 
ing  the  rotor  speed  and  overall  mass  flow. 

The  radial  momentum  equation  was  developed  from  the  Euler 
equations  by  assuming  that  entropy  and  enthalpy  do  not  change  along 
streamhnes  except  as  prescribed  across  the  blades  by  the  loss  and  devi¬ 
ation  correlations.  The  set  of  conservation  equations  are  simplified  by 
replacing  the  axial  momentum  and  the  energy  equation  with  AS  =  0  and 
AH  =  0.  The  remaining  two  equations,  continuity  and  radial  momen¬ 
tum,  are  solved  simultaneously  at  each  axial  station. 

Modifications  of  the  SLCC  Technique 

For  the  TEACC  simulation  to  be  responsive  to  the  local  change  in 
total  pressure  due  to  inlet  distortion  and  capable  of  modeling  transients, 
the  sLcC  must  be  restricted  to  a  small  axial  region  on  either  side  of  the 
blades.  Since  the  SLCC  is  a  subsonic  flow  solver,  it  needs  a  full  set  of 
boundary  conditions  around  its  domain.  The  SLCC  boundary  condi¬ 
tions  of  inlet  and  exit  curvature,  overall  mass  flow,  swirl  angle,  total 

inlet  Boundary  Conditions  Exit  Boundary  Condition 
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Fig.  4.  Modifications  to  SLCC  for  calculating  source  terms. 
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temperature,  and  total  pressure  are  calculated  from  the  TEACC  flow 
field.  Fig.  4  is  included  to  show  the  truncated  grid  in  which  the 
SLCC  is  restricted  to  operate. 

Although  the  boundary  conditions  are  of  the  same  type  in  this 
new  mode  of  operation  as  they  were  before  the  grid  was  truncated  to 
the  blades,  they  are  now  more  complicated  since  they  are  no  longer 
uniform  but  a  function  of  the  time  dependent  TEACC  flow  field 
around  the  blades. 

The  boundary  conditions  for  the  modified  SLCC  technique  are 
outlined  below  for  the  computational  domain  specified  in  Fig.  4. 
Inlet  and  exit  curvature  provided  by  the  3D  integrator  are  calculated 
from  a  circumferential  projection  of  the  flow  field  onto  a  circumfer¬ 
ential  slice  (axial-radial  plane).  This  is  effectively  the  same  as  cal¬ 
culating  streamline  curvature  from  only  the  axial  and  radial  velocity 
components  of  the  3D  transient  velocity  flow  field.  Curvature  is  cal¬ 
culated  as  a  function  of  local  velocity  and  acceleration.  Overall 
mass  flow  rate  as  a  boundary  condition  to  the  SLCC  is  calculated 
from  the  3D  flow  field  just  upstream  of  the  bladed  region  by  the 
integration  of  mass  for  each  circumferential  segment.  The  mass  flow 
rate  of  each  circumferential  segment  is  summed  together  to  calculate 
overall  mass  flow.  A  radial  distribution  of  swirl  angle  is  calculated 
from  the  TEACC  solution  at  the  blade  inlet  and  defined  as  the  arc-tan- 
gent  of  tangential  velocity  divided  by  axial  velocity.  The  total  tempera¬ 
ture  and  total  pressure  at  the  inlet  are  calculated  by  using  the  conserva¬ 
tion  variables,  ideal  gas  relations,  and  the  compressible  form  of  the 
stagnation  definitions  at  the  inlet  of  each  control  volume. 


Technique  to  Calculate  Turbomachlnerv  Source  Terms 

The  source  term  calculations  are  performed  after  the  SLCC  con¬ 
verges  to  a  steady-state  solution  through  the  bladed  region.  The  tech¬ 
nique  for  calculating  turbomachinery  source  terms  uses  control  vol¬ 
umes  within  the  bladed  region  and  applies  steady-state  conservation 
laws  across  each  control  volume.  Since  the  conservation  of  angular 
momentum  is  maintained  in  the  axisymmetric  solution  of  the  SLCC,  a 
radial  distribution  of  circumferential  velocity  vectors  are  produced. 
Cartesian  control  volumes  are  constructed  over  a  circumferential  seg¬ 
ment  of  the  bladed  region  from  streamlines  and  overall  blade  geometry 
with  velocities  and  pressures  known  on  all  surfaces  fix)m  streamline 
calculations.  A  simplifying  assumption  that  the  top  and  bottom  surfaces 
of  each  control  volume  are  streamsurfaces  is  incorporated  since  mass, 
momentum,  and  energy  does  not  cross  a  streamsurface.  The  tuiboma- 
chinery  forces  are  calculated  from  a  control  volume  analysis  and  the 
streamsurface  assumption  by  the  pressure  area  forces  and  the  inertial 
forces  (Fig.  5). 

METHODOLOGY  TO  IMPLEMENT  SOURCE  TERMS  INTO 
TEACC 

With  the  TEACC-supplied  boundary  conditions,  the  SLCC  pro¬ 
duces  an  axisymmetric  solution  through  the  bladed  region  in  the  form 
of  streamsurfaces  which  are  constructed  into  control  volumes  to  calcu¬ 
late  turbomachinery  source  terms.  A  methodology  for  distributing  the 
sources  radially,  circumferentially,  and  axially  through  the  grid  packed 
region  through  the  blade  was  developed. 


Radial  Interpolation  of  Turbomachlnerv  Source  Terms 

To  implement  these  source  terms  into  TEACC,  a  radial  interpola¬ 
tion  technique  was  derived.  A  radial  distribution  of  sources  is  con¬ 
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Fig.  5.  Turbomachinery  blade  force  calculations. 

structed  by  selecting  the  radius  at  the  center  of  each  SLCC  source  con¬ 
trol  volume  (Fig.  6a).  This  radius  is  nondimensionalized  by  the  radial 
extent  of  the  inner  and  outer  casing  in  the  axial  center  of  the  blade  (Fig. 
6b). 

A  spline  is  used  to  interpolate  these  source  terms  to  TEACC  along 
a  single  circumferential  segment  of  the  bladed  region.  A  radial  distribu¬ 
tion  of  TEACC  control  volumes  is  defined  through  the  bladed  region 
using  the  existing  grid  structure  as  shown  in  Fig.  6c.  The  radius  is  non¬ 
dimensionalized  by  the  radial  extent  of  the  inner  and  outer  casing  in  the 
middle  of  the  bladed  region.  A  radial  distribution  of  TEACC  source 
terms  is  acquired  by  interpolating  the  fixed  TEACC  volume  centers 
onto  the  function  of  sources  developed  earlier  ftom  the  SLCC  as  shown 
in  Fig.  6d. 

Circumferential  Distribution  of  Turbomachlnerv  Source 
Terms 

The  circumferential  distribution  of  turbomachinery  source  terms 
are  calculated  by  a  direct  application  of  the  SLCC  in  each  circumferen¬ 
tial  grid  segment.  For  a  distorted  inlet  upstream  of  the  compressor,  a 
circumferential  inlet  distortion  is  imposed  on  the  system,  as  in  Fig.  7a. 
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Fig.  6.  Interpolation  of  radial  distribution  of  source  terms. 
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Fig.  7.  Circumferential  distribution  of  source  terms. 


This  distortion  affects  the  adjacent  flow  field,  and  Fig.  7b  shows  that 
this  circumferential  distortion,  ignoring  viscous  effects,  will  persist  up 
to  the  compressor  inlet  (Longley,  1992).  The  compressor’s  perfor¬ 
mance  wiU  be  affected  by  this  distortion,  causing  a  nonuniform  circum¬ 
ferential  flow  field.  The  modeling  technique  for  acquiring  a  circumfer¬ 
ential  distribution  of  turbomachinery  source  terais  applies  the  SLCC 
separately  at  each  circumferential  segment,  as  shown  in  Fig.  7c.  The 
SLCC  is  sensitive  to  the  changing  demands  to  the  flow  field  because  it 
acquires  its  boundary  conditions  immediately  upstream  and  down¬ 
stream  of  each  circumferential  segment.  Therefore,  the  SLCC  interpo¬ 
lates  a  new  radial  distribution  of  turbomachinery  source  terms,  sensi¬ 
tive  to  the  changing  flow  field  for  each  circumferential  segment  as 
shown  in  Fig.  7d.  With  a  more  complex  inlet  distortion  pattern,  the 
number  of  circumferential  segments  would  be  increased  to  maintain  the 
high  fidelity  of  the  TEACC  simulation. 


over  the  top  40  percent  of  the  blade  while  a  double  circular-arc  con¬ 
struction  was  employed  for  the  bottom  60  percent  of  the  blade.  With  a 
hub-to-tip  ratio  at  the  rotor  inlet  of  0.5,  the  blade  sections  were  long 
enough  to  require  a  mid-span  damper  to  maintain  structural  integrity 
during  operation.  Rotor  IB  consisted  of  44  blades,  producing  a  moder¬ 
ate  solidity  of  1.3  at  the  rotor  tip. 

Grid  Development 

TEACC  requires  a  fixed  three-dimensional  grid  on  which  to 
resolve  the  conservation  equations.  A  grid  of  69  x  13  x  26  was  con¬ 
structed  to  model  Rotor  IB,  where  the  inlet  and  exit  of  the  blade  were 
defined  as  slanting  lines  in  Fig.  8.  The  grid  was  constmcted  with  verti¬ 
cal  grid  hnes  in  the  vicinity  of  the  screen  to  properly  model  the  experi¬ 
mental  distortion  screen.  Grid  lines  were  smoothly  packed  to  reduce 
numerical  losses  through  the  bladed  region  where  the  flow  was  known 
to  have  strong  gradients.  A  radial-circumferential  view  of  this  base  grid 
is  presented  in  Fig.  9,  where  a  cylindrical  right-handed  coordinate  sys¬ 
tem  has  been  used  with  uniform  circumferential  segments.  An  integer 
number  of  grid  segments  was  constructed  to  conveniently  model  a  90- 
deg,  one-per-revolution  inlet  distoition.  Each  segment  was  15  degs  in 
circumferential  extent  with  a  total  of  24  segments  within  a  circle  (6 
segments  in  each  quadrant). 

A  grid  density  investigation  was  conducted  by  doubling  the  grid 
independently  in  each  direction.  This  study  revealed  that  the  radial  dis¬ 
tribution  of  the  flow  field  was  unaffected  by  the  increase  in  grid  den¬ 
sity.  However,  the  error  in  overall  total  pressure  ratio  and  efficiency 
was  roughly  halved  by  doubling  the  grid  in  the  circumferential  direc¬ 
tion.  This  improvement  was  found  to  be  directly  related  to  halving  the 
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Axial  Distribution  of  Turbomachinerv  Source  Terms 

The  grid  was  packed  through  the  bladed  region  to  reduce  the 
numerical  error  from  strong  flow  gradients.  This  method  means  a  great 
deal  of  freedom  exists  in  how  the  sources  should  be  distributed  through 
the  bladed  region.  The  sources  are  distributed  conservatively  through 
the  bladed  region  by  requiring  that  the  sum  of  the  sources  distributed 
through  the  NPARC  grid  must  equal  the  sum  of  the  sources  developed 
in  the  SLCC.  A  simple  technique  to  investigate  this  problem  was  incor¬ 
porated  by  using  a  weighting  function  which  could  take  on  a  variety  of 
linear  shapes.  However,  a  uniform  weighting  function  was  found  to  be 
the  most  robust. 

VALIDATION  RESULTS 

Rotor  IB  (Seyler,  1967)  was  chosen  as  the  validation  vehicle 
because  it  offered  simplicity  in  the  number  of  stages  (1  blade  row),  and 
it  represented  a  compression  system  with  a  thorough  analysis  of  clean 
and  distorted  inlet.  Rotor  IB  is  a  high-performance  transonic  rotor  sim¬ 
ilar  to  those  found  in  modern  high-speed  aircraft.  The  rotor  was 
designed  with  a  multiple  circular-arc  blade  shape  which  was  applied 


Fig.  8.  Axial-radial  view  of  3D  grid  modeling  Rotor  IB  and 
test  facility. 
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Fig.  9.  Radlal-circumferentia!  view 
of  3D  grid  for  Rotor  IB. 


error  in  cross  sectional  area  and  mass  flow  by  doubling  the  circumfer¬ 
ential  grid. 

CLEAN  INLET  VALIDATION  OF  TEACC 

TEACC  was  tailored  to  simulate  Rotor  IB  by  calibrating  the 
streamline  curvature  code  (SLCC)  used  to  calculate  the  turbomachinery 
source  terms.  The  calibration  was  done  by  specifying  a  radial  distribu¬ 
tion  of  loss  and  exit  relative  flow  angle  from  experimental  data  as  a 
function  of  corrected  mass  flow  and  by  specifying  a  radial  distribution 
of  blockage  as  a  function  of  corrected  mass  flow. 

Calculated  overall  total  pressure  map  of  Rotor  IB  with  a  clean  inlet 
is  presented  in  Fig.  10  and  compared  with  experimental  data.  Four  cor¬ 
rected  speeds  (50,  70,  90,  and  100  percent)  were  available  for  compari¬ 
son,  with  symbols  presented  for  the  experimental  data  and  lines  for  the 
TEACC  solution.  In  an  overall  sense,  there  is  good  comparison 
between  TEACC  and  data.  The  maximum  percent  difference  in  total 
pressure  ratio  is  less  than  1 -percent  difference  for  each  characteristic. 

TEACC’ s  overall  efficiency  calculations  for  Rotor  IB  are  com¬ 
pared  with  experimental  data  in  Fig.  1 1.  The  overall  shape  and  charac¬ 
ter  compares  well  with  the  experimental  data.  Each  efficiency  curve 
has  a  peak  value  which  TEACC  fails  to  reach,  causing  a  maximum  per¬ 
cent  difference  (for  all  speeds)  between  TEACC  and  experimental  data 
of -2.9  percent  occurring  at  the  50  percent  corrected  speedline. 

Radial  Comparison  of  TEACC  with  Data  (Clean  Inlet) 

A  radial  comparison  was  made  with  exit  total  pressure  and  total 
pressure  ratio  as  presented  in  Figs.  12  and  13  for  the  100-percent  speed 
point  near  the  design  throttle  line.  The  overall  shape  is  quite  compK- 


Flg.  10,  Clean  Inlet,  overall  pressure  ratio  compared  with 
experimental  data. 
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Fig.  11.  Clean  inlet,  efficiency  compared  to  experi¬ 
mental  data. 


cated  for  total  pressure,  increasing  from  hub  to  tip  with  the  hub  calcula¬ 
tions  deviating  the  most  from  experimental  data  by  -1.3  percent  The 
tir*  and  overall  shane  are  in  nood  agreement  with  the  experimental  data.. 
The  total  pressure  ratio  is  compared  with  data  in  Fig.  13,  where  a  wide 
radial  ranne  in  total  pressnre  ratio  ean  he  seen.  There  is  nood  agreement 
(1.0  percent  difference)  between  TEACC  and  the  experimental  data. 
The  exit  st.atie  pressure  and  exit  IV/faeh  number  eonsistentiv  .support  the 
strong  distribution  in  total  pressure  ratio.  As  shown  in  Figs.  14  and  15, 
the  exit  Mach  number  and  the  exit  static  press.ure  a.re  in  nood  a.nree- 


100-percent  speedllne. 


Fig.  13.  Clean  Inlet,  total  temperature  ratio  on 
the  100-percent  speedllne. 


Fig.  14.  Clean  Inlet,  exit  Mach  number  on  the 
100-percent  speedllne. 


Fig.  15.  Clean  Inlet,  exit  static  pressure  on  the 
100-percent  speedllne. 


ment  with  the  experimental  data,  and  the  exit  Mach  number  is  fairly 
uniform,  but  the  exit  static  pressure  varied  consistent  with  the  exit  total 
pressure.  The  largest  differei^  for  both  exit  Mach  number  and  exit 
static  pressure  with  experimental  data  is  at  the  tip,  where  Mach  number 
is  underpredicted  by  2.7  percent  and  static  pressure  is  overpredicted  by 
1.9  percent.  The  SLCC  operate^n  the  relative  reference  frame  where 
the  relative  Mach  number  was  observed  to  be  transonic.  The  TEACC 
simulation  remains  in  the  absolute  reference  frame  where  Mach  num¬ 
ber  remained  subsonic  (Hg.  14). 


Experimental  Data 


DISTORTED  INFLOW  VALIDATION  OF  TEACC 

Rotor  IB  was  tested  with  an  inlet  screen  to  quantify  the  effects  of 
inlet  flow  distortion  on  its  performance.  Because  the  distortion  data 
provided  with  Rotor  IB  was  based  on  a  90-degree,  one-per-revolution 
screen,  a  simulation  of  the  screen  was  used  in  the  TEACC  simulation. 
For  this  study,  TEACC  was  compared  at  three  different  conected 
speedlines  for  a  distorted  inlet  A  single  distortion  point  on  the  100-per- 
cent  speedline  was  investigated  in  great  detail  because  radial  and  cir¬ 
cumferential  data  were  available.  Overall  performance  of  the  distorted 
compressor  is  compared  to  experimental  data  as  well  as  radial  and  cir¬ 
cumferential  distributions. 

Distortion  Screen  Modeling 

To  model  the  distortion  screen,  the  porous  wall  boundary  condition 
(Cooper,  K.,  Jones,  R.,  and  Kincade,  B.  “NPARC  Porous  Wiill  Bound¬ 
ary  Conditions.”  private  communication  and  notes.  May  1991;  Pinker, 
1967)  was  used.  The  boundary  condition  is  semi-empirical  in  that 
experimental  data  has  been  taken  across  a  variety  of  screens  to  estab¬ 
lish  a  loss  in  total  pressure  as  a  ftinction  of  Mach  number,  porosity,  and 
Reynolds  number.  The  value  of  total  pressure  at  the  exit  of  the  screen 
was  calculated  from  the  continuity,  energy  equations  (total  enthalpy 
equal  to  a  constant),  and  an  empirical  pressure  loss  coefficient.  This 
general  technique  to  simulate  the  effect  of  a  distortion  screen  was 
incorporated  for  aU  inlet  distortion  calculations. 

A  single  screen  with  a  90-deg,  one-per-revolution  pattern  was 
designed  to  give  a  classical  circumferential  inlet  total  pressure  inlet  dis¬ 
tortion  to  Rotor  IB.  At  the  100-percent  speed,  the  experimental  data 
showed  that  the  screen  produced  a  total  pressure  loss  of  15  percent 
from  the  clean  inlet  which  was  duplicated  for  the  TEACC  simulation. 
Investigations  were  conducted  through  the  middle  of  the  four  circum¬ 
ferential  quadrants  where  the  experimental  data  was  taken.  The  screen 
was  located  in  quadrant  A  (Fig.  16).  The  compressor  rotates  clockwise, 
causing  air  to  swirl  in  the  direction  of  increasing  quadrant  letter 
(A,  B,  C,  and  D). 

Illustrated  in  Fig.  17  are  the  radial  characteristics  of  the  distor¬ 
tion  screen  simulation  and  comparison  to  experimental  results 
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Fig.  17.  Inlet  total  pressure  distortion  produced  by  the 
screen  simulation  at  100-percent  corrected  speed. 

with  experimental  data  in  Fig.  17a  and  TEACC  results  in  Fig.  17b.  The 
character  of  the  flow  field  behind  the  screen  has  certainly  been  cap¬ 
tured,  with  Quadrant  A  (90-deg  distortion)  clearly  separated  from  the 
other  three  quadrants.  The  individual  pressure  drops  for  Quadrants  B-D 
calculated  by  the  screen  simulator  are  clustered  together  and  represent 
the  mean  of  the  data  in  these  quadrants.  Not  enough  experimental  data 
was  available  to  adequately  explain  the  scatter  in  the  experimental  data 
for  Quadrants  B-D.  A  possible  explanation,  however,  is  that  the  scatter 
is  a  result  of  unsteadiness  in  the  flow  field  upstream  of  the  distortion 
screen. 

Overall  Rotor  IB  Performance  Comparison  (Distorted  Inlet) 

An  overall  total  pressure  ratio  performance  map  for  Rotor  IB  is 
presented  in  Fig.  18  for  three  corrected  speeds  of  50,  70,  and  100  per¬ 
cent  The  clean  speedlines  are  presented  to  give  a  proper  orientation  of 
the  distorted  data,  with  stall  depicted  for  each  speedhne  at  its  lowest 
corrected  mass  flow  rate.  A  conservative,  but  simple,  stalling  criterion 
was  imposed  on  the  TEACC  simulation.  When  any  one  of  the  circum¬ 
ferential  segments  acquired  a  corrected  mass  flow  which  exceeded  stall 
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Fig.  18.  Distorted  Inlet,  overall  performance  map  with  total  pressure 
ratio  at  100-percent  corrected  speed. 


Fig.  16.  Screen  location  and  labeling. 


for  the  clean  inlet,  the  compressor  was  considered  stalled.  Likewise,  the 
TEACC  simulation  was  halted  when  any  one  of  the  circumfCTential 
seements  exceeded  the  clean  inlet  data  on  the  choked  end  of  the  speed¬ 
lines,  These  two  restrictions  define  the  distortion  calculation  limits  pre¬ 
sented  for  each  speed  in  He.  18, 

For  the  two  low-speed  cases,  the  clean  and  distorted  data  lay  close 
together.  The  extent  in  corrected  mas.s  flow  rate  is  about  the  same  for 
7the  simulation  and  the  experimental  distortion  data.  The  maximum  per- 
^  cent  difference  between  TEACC  results  and  experimental  data  at  the 
lower  speeds  was  approximately  2  percent  The  lOO-percent  speed  is 
more  interestine  because  the  experimental  data  presented  a  region 
where  the  compressor  intermittently  stalled  or  remained  stable  with  the 
imposed  inlet  distortion  screen.  The  TEACC  simulation  predicted  the 
corrected  mass  flow  rate  differently  from  reported  experimental  data  by 
a  maximum  of  3  percent.  The  TEACC  simulation  maps  out  a  ranee  of 
corrected  mass  flow  which  includes  approximately  the  middle  of  the 
intermittent  stall  reeion. 

Overall  efficiency  is  compared  to  experimental  data  in  Fig.  19.  In 
all  cases  TEACC  predicts  an  efficiency  that  is  higher  than  reported 
experimental  data.  The  maximum  percent  difference  between  TEACC 
and  the  experimental  data  is  located  near  the  stall  or  choked  portion  of 
each  corrected  speed.  The  maximum  percent  difference  between 
TEACC  results  and  data  are  6.9, 4.5,  and  5,9  percent  difference  for  50-, 
70-,  and  1 00-percent  speeds,  respectively. 

Radial  Comparison  of  Compressor  with  Experimental  Data 
^Distorted  Inlet) 

This  investigation  now  turns  to  an  examination  of  TEACC  results 
compared  to  experimental  data  behind  the  compressor.  Exit  total  pres¬ 
sure  (defined  by  the  ratio  of  exit  total  pressure  over  the  average  screen 
inlet  total  pressure)  versus  compressor  radius  is  plotted  with  experi¬ 
mental  data  in  Fig.  20a  and  TEACC  results  in  Fig.  20b.  Pressure  ratio  is 
characterized  by  a  division  between  Quadrant  A  and  the  other  three 
quadrants.  TEACC  does  a  good  job  of  matching  the  character  of  the 
exit  total  pressure  in  Quadrant  A  with  experimental  data.  TEACC 

Experimental  Data 


results  agree  with  the  experimental  data  in  the  general  character  of  the 
other  three  quadrants  and  identifies  Quadrant  B  with  the  largest  exit 
total  pressure.  The  experimental  data  also  confirms  the  TEACC  results 
and  shows  that  exit  total  pressure  decreases  steadily  from  Quadrant  B 
to  Quadrant  D.  Although  the  shapes  of  the  other  three  quadrants  are 
moving  in  the  correct  direction,  the  extent  and  scatter  with  pressure 
ratio  is  mis.sing.  The  TEACC  results  in  Quadrant  A  compare  well  to  the 
experimental  data,  except  for  the  hub  point,  which  causes  the  maxi¬ 
mum  percent  difference  with  pressure  ratio  equal  to  6.15  percent. 
TEACC  results  compared  to  experimental  data  are  too  low  for  Quad¬ 
rant  B  (6.82  percent)  and  too  high  for  Quadrants  C  (5.86  percent)  and  D 
(8.64  percent).  The  highest  maximum  percent  difference  is  within  8.7 
percent  in  Quadrant  D. 

The  final  comparison  with  experimental  data  is  also  behind  the 
compressor  (same  reference  location  as  exit  total  pressure).  The  ratio  of 
total  temperature  at  the  compressor  exit  divided  by  the  average  total 
temperature  at  the  screen  inlet  is  compared  to  experimental  data  with 
data  presented  in  Fig.  21a  and  TEACC  results  in  Fig.  21b.  Quadrant  A 
is  again  characterized  as  having  a  greater  slope  than  the  other  three 
quadrants,  and  TEACC  does  a  good  job  identifying  this  overall  flow 
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Fig.  19.  Distorted  inlet,  overall  performance  map 
with  efficiency. 
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Fig.  20.  Distorted  Inlet,  exit  total  pressure  on  100-percent  speedline. 
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Fig.  21.  Distorted  Inlet,  exit  total  temperature  on  100-percent  speedline. 


character.  TEACC  correctly  identifies  that,  in  the  nondistorted  region, 
the  exit  total  temperature  was  the  highest  in  Quadrant  B,  with  the  other 
quadrants  decreasing  in  value  from  Quadrants  B-D.  TEACC’ s  calcula¬ 
tions  for  the  total  temperature  at  the  compressor  exit  have  a  greater 
spread  in  Quadrants  B-D  than  for  the  exit  total  pressure.  Comparisons 
also  show  that  the  hub  of  Quadrant  A  is  not  as  steep  a  slope  as  the 
experimental  data,  but  the  tip  matches  the  data  very  well.  Quadrants  C 
and  D  match  the  data  well,  within  1.81  percent  and  2.19  percent, 
respectively,  and  even  Quadrant  B  matches  the  data  weU,  except  at  a 
single  tip  point  The  largest  percent  difference  of  TEACC  with  experi¬ 
mental  data  for  exit  total  temperature  is  6.60  percent. 


EFFECT  OF  SWIRL  ON  THE  COMPRESSOR  PERFOR¬ 
MANCE  WITH  INLET  DISTORTION 

The  circumferential  velocity  within  the  flow  field  caused  by  the 
high  rotor  speed  compressor  blade  is  a  direct  result  of  the  compressor 
blade  doing  work  on  the  flow.  The  blade  causes  a  strong  circumferen¬ 
tial  velocity  which  establishes  the  location  of  the  largest  static  and  total 
pressure  regions.  The  compressor  exit  total  pressure  is  presented  in  Fig. 
22.  Notice  that  the  highest  value  in  total  pressure  is  not  in  the  center  of 
the  distorted  region,  but  has  shifted  clockwise  and  has  its  largest  value 
near  the  tip.  Displayed  in  Fig.  23  is  exit  total  temperature,  where  it  has 
shifted  clockwise  but  is  not  as  pronounced  as  the  exit  total  pressure. 
Figure  24  presents  the  static  pressure  which  resembles  that  of  the  previ¬ 
ously  presented  total  pressure.  The  highest  static  pressure  is  close  to  the 
interface  of  Quadrant  A  and  B.  The  velocity  vectors  can  be  converted 
to  streamlines  which  show  the  strong  effect  of  swirl  downstream  of  the 
compressor. 


Fig.  22.  Exit  total  pressure  distribution 
produced  by  90-deg,  one-per- 
revolutlon  screen  on  the  100- 
percent  speedllne. 


Fig.  23.  Exit  total  temperature  distribution 
produced  by  90-deg,  one-per- 
revolutlon  screen  on  the  100- 
percent  speedllne. 


Fig.  24.  Exit  static  pressure  distribution 
produced  by  90-deg,  one-per- 
revolutlon  screen  on  the  100- 
percent  speedllne. 


The  streamlines  are  shown  in  Fig.  25.  A  schematic  of  the  distortion 
screen  is  presented  to  indicate  inlet  distortion  (Fig.  25a);  Quadrant  A  is 
defined  in  the  middle  of  the  distorted  region  and  Segm^t  D  is  chosen 
to  represent  a  segment  just  before  the  rotor  begins  to  sweep  through  the 
distortion  in  Quadrant  D.  An  important  distinction  is  that  Quadrant  A  is 
in  the  distorted  region,  and  Segment  D  is  in  the  undistorted  region.  The 
upper  right-hand  picture  in  Fig.  25b  shows  in  simple  terms  how  Quad¬ 
rant  A  and  Segment  D  behave  on  the  compressor  map.  First,  Segment 
D  operates  in  the  clean  inlet  region  of  the  compressor,  producing  a 
pressure  ratio  with  a  high  corrected  mass  flow,  but  Quadrant  A  operates 
through  the  distorted  screen  where  the  total  pressure  in  front  of  the 
compressor  is  low,  causing  a  reduction  in  corrected  mass  flow.  Quad¬ 
rant  A  operates  higher  on  the  compressor  map  closer  to  stall,  which 
means  that  approximately  75  percent  of  the  compressor  operates  near 
choke  while  25  percent  of  the  compressor  operates  near  stall.  The  por¬ 
tion  of  the  compressor  operating  behind  the  distortion  produces  more 
work  than  any  other  part  of  the  compressor. 

The  streamlines  in  the  bottom  view  of  Fig.  25c  tell  more  about  the 
compressor  operation.  Segment  D  describes  the  streamlines  to  the  right 
of  the  distorted  region  while  the  blade  turns  the  flow  strongly  toward 
and  behind  the  distorted  region.  The  streamlines  through  the  middle  of 
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c.  Streamlines  through  Rotor  IB 
Fig.  25.  Streamlines  showing  pressure  convected  with  the  flow 
on  the  100-percent  speedllne  with  Inlet  distortion. 


the  distorted  region  are  represented  by  Quadrant  A.  Again,  the  com¬ 
pressor  strongly  swirls  the  flow,  but  the  total  pressure  increase  through 
the  distorted  region  is  greater  than  that  of  Segment  -D.  The  flow  field  at 
the  compressor  exit  is  convected  clockwise  by  the  swirling  flow  to  an 
adjacent  part  of  the  compressor,  and  the  streamlines  in  Segment  P 
show  that  the  lower  total  pressure  fluid  is  convected  behind  the  com¬ 
pressor  in  the  distorted  region.  This  is  a  further  indication  of  the  clock¬ 
wise  transport  of  total  pressure  attributable  to  the  swirling  flow,  and  an 
indication  that  the  total  pressure  behind  the  compressor  is  not  uniform, 
but  increasing  in  the  direction  of  flow  rotation.  The  performance 
increase  which  begins  at  the  distorted  region  increases  the  static  pres¬ 
sure  behind  the  compressor  and  convects  total  temperature  and  total 
pressure  with  the  flow  to  an  adjacent  point  behind  the  compressor.  This 
action  causes  the  adjacent  point  to  further  backpressure  the  compressor 
in  this  region,  subsec^uentlv  increasing  its  performance  and  moving  that 
region  closer  to  stall.  This  increase  implies  that  the  swirl  will  cause  the 
stronge^st  backpressuring  of  the  distorted  region  to  be  transported  clock¬ 
wise  behind  the  adjacent  non-distorted  region  (intersection  of  Quadrant 
A.  and  B),  increasing  perfoimance  in  this  region  and  producing  the  larg¬ 
est  total  pressure  found  anywhere  in  the  compressor,  as  shown  in  Fig. 
22.  Therefore,  the  total  pressure  increases  clockwise  through  the  dis¬ 
torted  region. 

The  transport  of  total  pressure  and  total  temperature  clockwise 
through  the  compressor  can  be  viewed  in  a  slightly  different  way 
Fig.  26,  Here,  a  radial  mass-wsdghted  average  of  total  pressure 
and  total  temperature  is  presented  versus  the  circumferential  extent  at 
the  upstream  and  downstream,  reference  locations.  Inlet  total  pressure 
(Fig.  26b)  to  the  compressor  can  be  observed  to  have  two  different  uni¬ 
form.  regions,  one  behind  the  distortion  screen  (Quadrant  A)  and  a 
higher  value  flnough  the  other  three  quadrants.  The  compressor  exit 
total  pressure  incre^ases  ste^adily  through  the  distorted  region.  y,dth  its 
highest  value  at  the  beginning  of  Quadrant  B.  The  ratio  of  total  pressure 
(Fig.  26a)  also  increases  steadily  in  a  clockwise  direction  through  the 
distorted  region,  due  to  the  compressor  transportation  of  total  pressure 
clockwise  through  the  distorted  region.  This  increases  the  .static  pres¬ 
sure  behind  an  adjacent  portion  of  the  compressor  in  the  clockwise 


direction,  backpressuring  this  portion  of  the  compressor  and  producing 
a  performance  increase.  This  increased  performance  shows  itself  as  an 
increase  in  total  pressure  ratio. 

A  similar  story  can  be  constructed  with  total  temperature  at  the 
compressor  inlet  and  exit.  The  inlet  total  temperature  to  the  compressor 
remains  uniform  in  the  presence  of  the  screen,  as  shown  in  Fig,  26d. 
The  exit  total  temperature,  like  the  exit  total  pressure,  steadily  increases 
clockwise  through  the  distorted  region.  This  means  the  total  tempera¬ 
ture  ratio  in  Fig.  26c  also  increases  through  the  distorted  region.  Total 
temperature  ratio,  just  like  total  pressure  ratio,  reaches  its  largest  value 
at  the  interface  of  Quadrant  A  and  Quadrant  B  because  of  the  swirl 
effect  through  the  distorted  region. 

SUMMARY/CONCLUSIONS 

TEACC  has  been  developed  to  model  inlet  distortion  through  the 
use  of  semi-actuator  disk  theory.  This  technique  allows  efficient  com- 
putationalgridsthatprovideapracticalsolutionforinletdistortionproblems. 
TEACC  was  constructed  by  developing  a  methodology  to  merge  two 
proven  technologies,  NPARC  and  a  SLCC,  into  a  single  turbomachin¬ 
ery  simulation. 

TEACC’ s  predictions  with  experimental  data  (Rotor  IB)  for  a 
clean  inlet  compared  within  1  percent  for  total  pressure  ratio  and  within 
3  percent  difference  for  efficiency.  A  radial  comparison  of  TEACC 
with  experimental  data  also  produced  good  results. 

The  distortion  was  produced  with  a  90-deg,  one-per-revolution 
screen  simulated  by  a  boundary  condition  provided  by  NPARC. 
TEACC  simulated  Rotor  IB  over  three  speedlines  and  produced  a  total 
pressure  ratio  which  compared  within  3  percent  of  experimental  data 
and  an  efficiency  which  compared  within  7  percent  difference  of  exper¬ 
imental  data. 

A  radial  comparison  of  TEACC  results  to  experimental  data 
revealed  that  TEACC  models  the  overall  character  of  the  compressor 
well  and  does  a  particularly  good  job  of  predicting  the  magnitude  and 
shape  of  exit  total  temperature  and  exit  total  pressure  in  the  distorted 
region  (Quadrant  A). 

A  detailed  investigation  of  the  compressor  flow  field  was  con¬ 
ducted  to  evaluate  the  effects  of  the  distortion  screen  simu- 
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lation  and  to  compare  these  results  with  experimental  data. 
The  study  revealed  that  strong  swirl  convects  total  tempera¬ 
ture  and  total  pressure  with  the  flow  in  the  direction  of  rotor 
rotation.  This  increase  in  total  pressure  produced  an 
increase  in  static  pressure  in  an  adjacent  part  of  the  com¬ 
pressor  which  backptessured  the  compressor,  producing  an 
increase  in  performance.  Therefore,  the  greatest  total  pres¬ 
sure  is  produced  in  transition  from  distorted  to  undistorted 
flow  regions. 
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The  TEACC  development  effort  to  simulate  Rotor  IB 
with  inlet  distortion  has  been  very  successful.  The  TEACC 
technique  was  developed  for  a  single  Rotor  and  expected  to 
be  immediately  expandable  to  multiple  blade  rows.  TEACC 
is  immediately  recognized  as  a  time  dependent  solution 
technique  which  can  be  extended  to  investigate  surge  and 
rotating  stall  instabilities.  TEACC  is  an  ideal  environment 
to  investigate  inlet  con^atibility  issues  with  the  inlet  and 
the  compressor  tightly  coupled. 
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